Complementarity Reveals Bound Entanglement of Two Twisted Photons 
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We witness for the first time the generation of bound entanglement of two photon qutrits, whose 
existence has been predicted by the Horodecki family in 1998. Detection of these heavily mixed 
entangled states, from which no pure entanglement can be distilled, is possible using a key concept 
of Nature: complementarity This captures one of the most counterintuitive differences between a 
classical and quantum world, for instance, the well-known wave-particle duality is just an example 
of complementary observables. Our protocol uses maximum complementarity between observables: 
the knowledge about the result of one of them precludes any knowledge about the result of the 
other. It enables ample detection of entanglement in arbitrary high-dimensional systems, including 
the most challenging case, the detection of bound entanglement. For this we manipulate "twisted" 
twin photons in their orbital angular momentum degrees of freedom. Our experimentally demon- 
strated maximum complementarity protocoF is very general and applies to all dimensions and 
arbitrary number of particles, thus enables simple entanglement testing in high-dimensional quan- 
tum information and opens up the quest of understanding the meaning of this type of entanglement 
in Nature. 



The boundary between a classical and quantum world is 
the domain of mixed quantum states. Quantum entan- 
glement of mixed states is much more complex than that 
of pure states and subject of intense research. In 1998 
it was found that entanglement of mixed states is not 
"fiat", but has an intriguing structure!!' Entangled quan- 
tum states can be classified into two distinct types, those 
that can and those that cannot be distilled into pure 
entangled states using statistical local operations and 
classical communication. Undistillable entangled states 
are called bound entangled. Bound entanglement can 
occur in bipartite systems with dimensions larger than 
(i = 2 or in the case of more particles. These two dif- 
ferent types of bound entanglement differ considerably. 
The case of more than 2 entangled qubits has recently 
been discovered experimentally in photonic multipartite 
systems j^M^ trapped ions NMRI^ further, there IS one 
study in a continuous- variable context.!^ These results 
on multipartite qubit bound entanglement became al- 
ready interesting for certain quantum information tasks 
such as steering^ and reduction of the communication 
complexity We investigate here the case of bound entan- 
glement of only two photonic qutrits {d = 3) using the or- 
bital angular momentum degree of freedom of light; this 
is the simplest case of bound entanglement and complica- 
tions such as occurring for multipartite system do n ot 
occur. The orbital angular momentum of photon is 
used in a number of studies to obtain high-dimensional 
bipartite entangled quditsj^^H^^ apart from fundamental 
questions such as those discussed here, the system shows 
promise in quantum cryptograph}/^^ due to improved se- 
curity, increased resistance against noise, and higher bit 
rateJI^H^ compared to qubits. In particular for such 



higher-dimensional qudit systems, it is well known that 
it is extremely hard to check whether a given state, es- 
pecially if it is mixed - the daily situation in a labora- 
tory - is separable or entangled. There are several oper- 
ational criteria available to detect entanglement, where 
the most famous and powerful one is the Peres-Horodecki 
criterion, which is based on the partial transposition in 
one subsystem .^^^^^ After partial transposition, either all 
eigenvalues remain positive (PPT), or in the other case, 
if at least one eigenvalue is found to be negative, the state 
is entangled. The transposition of a state is synonymous 
to time reversal, so PPT tests if a state where time's 
arrow is reversed for one of its partitions is still a phys- 
ical state: obviously, if the state is separable, this must 
be the case. For two qubits, the reverse is also true and 
PPT is a sufficient criterion for separability. This argu- 
ment fails already for two qutrits, where it turns out that 
bound entangled qutrits are PPT. This brings us back to 
distillability: it was found that it is actually the positiv- 
ity under partial transposition, which makes bound en- 
tangled states undistillable .'Sin the experiment described 
here, in a first step, we produce photon states that are 
bound entangled in their orbital angular momentum de- 
gree of freedom and verify the positivity under partial 
transposition via state tomography. In a second step we 
apply the maximum complementarity protocol using sets 
of observables, for which their eigenbases are mutually 
unbiased, to witness directly the inseparability of two 
bound entangled qutrits. 

The Maximum Complementarity Protocol: Con- 
sider the following scenario of a source producing two- 
qudit states p G C^^^, namely quantum states with d 
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FIG. 1. This figure illustrates the steps involved in the com- 
plementarity protocol for entanglement detection of bipartite 
qudits: Alice and Bob determine the experimental probabil- 
ity (coincidence counts) for each detection state, and that for 
all (mutually unbiased) bases. It is then simply a matter of 
summing up these probabilities, if the sum is > 2, the state 
is entangled. Note that both Alice and Bob are allowed to 
relabel their measurement outcomes per basis, in order to op- 
timize the detection ability. 

degrees of freedom per qudit. Both experimenters, Al- 
ice and Bob, can choose among k different observables. 
What is the best strategy for Alice and Bob to detect 
the inseparability? The most striking difference between 
entanglement and separability are revealed by correla- 
tions in different basis choices. A fully correlated sys- 
tem is a physical system for which we can predict with 
certainty the outcome of a second measurement when 
we know the outcome of a first measurement, opposite 
to the other extreme case when knowledge of the first 
measurement outcome does not reveal any information 
about the second measurement outcome. This we would 
call a fully uncorrelated system. Let us quantify this 
statement via a correlation function for two observables 
given in the spectral decomposition, and mea- 
sured by Alice and Bob (with respective eigenvectors 
{\ik)} = {|0/c), ...,\d- Ik)}) by summing all joint prob- 
abilities PAk,Bk{h^h) when both parties have the same 
outcome ik 

d-l d-1 

Ca^^b^ = ^PAk,BA^k^ik) = '^Tr{\ikik){ikik\ Pd) 

(1) 

Here, we allow Alice and Bob to relabel their outcomes 
such that C gets maximal. Indeed the above function is a 
mutual predictability, since if the state is fully correlated 
then the outcomes can always be labelled such that C = 1 
and if fully uncorrelated then any relabeling can only give 
C = 1/d since all d outcomes are equally likely. However, 
this function does not tell us anything about entangle- 
ment, since clearly the socks of Prof. Bertlmann,^^ who 
has the habit of wearing differently coloured socks, are 
also fully correlated. Therefore, in a second step, Alice 
and Bob will now use the fundamental concept of com- 



plementarity and choose a second set of observables that 
are mutually unbiased to the first choices of observables 
Ai and Bi, namely the observables A2 and B2. One way 
to phrase Bohr's complementarity of two observables Ai, 
A2 is to say that they are non-degenerate; i.e., they do 
not share any common eigenvector. This implies that the 
uncertainties of the outcomes are bounded by the scalar 
product of the eigenvectors of both observables and it is 
maximal if (and only if) all scalar products satisfy 

\{in\jm)\^ = 2 Vi,jG{0,l,...,rf-l}. (2) 

Thus, predictabilities of a product state |OiOi) are given 
by 

-PAi,Bi(il,«l) = <5i,0 

PA2,B,ii2,i2) = Tr(|i2i2)(i2i2|0i0i)(0i0i|) 

^'>£\\(^^.\(^^=}., (3) 

1 1 

d d 

which holds for all i and we obtain for the correlation 
functions Cai,Bi = 1 and Ca2,B2 = ^- For a third choice 
of mutually complementary observables we obtain again 
^. In general, by adding up m correlation functions we 
obtain 

m ^ 

Im:=J2^A,,B,<l + {m-l)- (4) 

k=l 

that has to hold for all pure separable states. In Ref.l^we 
have shown that it is also valid for all mixed separable 
states. Therefore, 1^ serves as a detection function of 
entanglement if and only if the bound is not satisfied. 
The existence of complete sets of mutually unbiased bases 
in arbitrary dimensions is an open problem, but in the 
case of prime-power dimensions d it is known to be (i + 1 , 
leading to 

_ d+i 

I,^^^'^J:CA.^B.^2. (5) 

k=l 

It was shown that this detection criterion is very pow- 
erful as it e.g. detects all the entanglement of isotropic 
states and more entanglement of multiparticle Aharanov 
state than the powerful criteria of Refs. pMMlj^ Let us 
remark that our method of maximum complementarity 
provides a test for non-separability of a state whereas 
Bell inequalities can serve as an entanglement test but 
are designed to test for local realism. Indeed, so far only 
multipartite bound entangled states were found to vio- 
late a Bell inequality while the question is still open for 
bipartite bound entangled states. 

Now, we demonstrate both theoretically and experimen- 
tally that is also capable of detecting bound en- 
tanglement for a certain class of so-called magic simplex 
states or Bell-diagonal states in 3 x 3 dimensions. 
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FIG. 2. Experimental setup for generation of bound- 
entangled bipartite qutrits. Photon pairs are created by 
downconversion (type-I) of 413 nm UV photons in a PPKTP 
crystal, momentum- filtered by pinhole (PH), and split proba- 
bilistically using a beamsplitter (BS). Rotation (ROT) in the 
orbital angular momentum superposition mode-space for pro- 
duction of the Bell states Pk,i (gray box) via W^fc^/s,z^/^ is 
performed using a spatial light modulator (SLM). Projective 
measurements (PROJ) are done by appropriate mode trans- 
formation and subsequent imaging onto the core of a single 
mode fibre (SMF) as indicated by Ha/b- We have imple- 
mented both operations on the same SLM. The photons are 
guided to an avalanche photo diode (APD), and detection 
events belonging to a single photon pair are post-selected by 
coincidence detection. 



For the experimental test, we have chosen to use 
orbital-angular- momentum^ ^ (GAM) entangled photons 
generated by spontaneous parametric downconversion 
(SPDC). GAM entanglement^^ is one implemen- 

tation of spatial entanglement which is particularly in- 
tuitive because the quantum correlations of the photon 
pair are determined only by conservation of orbital angu- 
lar momentum during pair generation: The GAM £a and 
is of the downconverted photons must fulfil + ^ b=0 
if the pump is a flat (Gaussian) beam with £ = 0.'^^'^ 
In the GAM basis the two-photon qutrit state as 

produced by the crystal is 
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{|-1,+1) + |0,0) + | + 1,-1)}. (6) 



In the following, we denote the states {\i = — 1), |^ = 
0),\£ = +1)} by {|0), |2)}. As a starting point 
we take the maximally entangled Bell state Po,o = 
\'^SPDc){'^SPDc\ and, via applying the unitary Weyl 
operators Wk^i := Xln^o ^ ^ \n){n + /| in one sub- 
system (e.g. on the photon of Alice), we can synthe- 
size the (P — 1 maximally entangled Bell states Pk^i = 
Wk,i 1 Po,o Wl I <S) t. Any convex combination of these 
(P Bell states forms a so called ''magic simplex^'^yV := 
{pd = Ylk~i=o^k,iPk,i\^k,i > ^,Y.f^Lo^k,i = !}• This 
reduced state space of locally maximally mixed states 
admits a simple geometrical representation (the vertices 
are the cP Bell states) and has been showri^^M^ to be 



powerful in addressing inseparability issues and quantum 
information theoretic questions. 

For our purpose let us reduce the dP — 1 dimensional pa- 
rameter space to three (four in case of d > 3) parameters 
Qi and consider the following family of states: 
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This family also includes for d = 3 the one-parameter 
Horodecki-state, the first found bound entangled state. ^1 
Namely, for qi ^2 = -|f, ^3 = with 

A G [0, 5]. This state is PPT for A G [1, 4] and was shown 
to be bound entangled for A G {3,4]. 

Generation of the states {d = 3): By expressing the 
totally mixed state ILg as the sum of all Bell states ^ Pk^i 
we find that all 9 Bell states need to be mixed to synthe- 
size the state ps. Note that there are 72 unitary equiva- 
lent possibilities^ to generate ps that we will exploit to 
deduce the error (see supplementary information). We 
introduce a measurement-based scheme to pro duce th ese 
mixed states similar to the authors in Refs.l^M^ ^^ I ^^ * for 
polarization but in our case for GAM qutrits using spa- 
tial light modulation, see Fig. |2] As described above we 
have to apply single-qutrit rotations (e.g., on qutrit A) 
using the Weyl matrix Wk^i to transform Po,o into any of 
the 8 other maximally entangled Bell states Pj^^i. This 
operation is implemented on the spatial light modula- 
tors (SLMs, see Fig. [2|, which allows us to generate the 
mixed state by time-multiplexing of the rotation opera- 
tors Wk,i for a particular choice of q^. Gur method is 
physically equivalent to tracing over a separate degree of 
freedom to create decoherence, such as the spectral one, 
a method which has also been used to generate mixed 
state^^^^ or stochastic rotation with an additional opti- 
cal element . ^^ l ^^l^^ l In our case we extend this to obtain 
full control over the high-dimensional mixture by using 
retroactive mixing: we record photon counts for each of 
the states Pk^i, and form the incoherent summation in a 
computer afterwards. This allows us to explore the magic 
simplex W completely, while still obtaining exactly the 
same result as if mixing would happen during the photon 
counting time. Projective measurements are performed 
by mode-conversion on the spatial light modulator and 
imaging onto a single- mode fibre. 

State Tomography and Magic Simplex: We first 
perform quantum state tomograph}/ of the bound en- 
tangled state. We keep qs = —0.5776 fixed, this allows 
us to visualize quantum state properties in 2D: Fig. [3] 
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FIG. 3. Experimental (a) and theoretical (b) 2D slice 

92} through the magic simplex for fixed g3 = -0.5776. All piG. 4. Experimental (a) and theoretical (b) real parts 

coloured points correspond to states having positive semidef- of the density matrices of the bound entangled state {qi = 

inite eigenvalues, therefore representing physical states. The -0.07, ^2 = -1.73, ^3 = -0.5774). The fidelity of the experi- 

blue area (curved region) covers the range of states with a pos- mental state is 0.98, the complementarity protocol applied to 

itive partial transpose (PPT). The red triangular area indi- the experimental state gives 2 - /d+i = -0.032 ± 0.003 and 

cates where the maximum complementarity protocol applied the minimum eigenvalue of the partially transposed state is 

to the experimentally generated states (a) or theoretical states Min [eig(p^^)] = +0.0101 ± 0.0017; the state is bound en- 

(b) is greater than 2, thus detecting entanglement. States tangled. It is expected to be purely real; the experimental 

where all three conditions are fulfilled are bound entangled. imaginary parts are all below 10%. The axes labels indicate 

We see that the experimental and theoretical geometries agree the 0AM quantum numbers i. 
very well. 



shows experimental (a) and theoretical case (b), where 
each coloured pixel {^1,^2} corresponds to a state, blue 
if PPT (curved region), and red if the state violates the 
maximum complementarity protocol, 2 — < (tri- 
angle). In regions where all conditions are fulfilled, we 
have non-separable states with a positive partial transpo- 
sition, i.e. bound entanglement. We find that, compared 
to the total area of physical states, the bound entan- 
gled states occupy a significant space in the chosen slice 
through the magic simplex; the extension along turns 
out to be even much larger, i.e., the bound entangled 
states occupies a rather large volume. Since the regions 
in Fig. [3] are continuous ("they do not have holes"), it is 
obvious from the geometry that states within the region 
where all criteria apply are bound entangled. However, 
it is important to check that the conditions are also ful- 
filled significantly.^^ In contrast to the polarization qubit 
case,^ in spatial entanglement, the most important er- 
rors are due to unwanted rotations of the state due to 
wavefront aberrations. To quantify these errors, we have 



determined for all 72 unitary-equivalent states of ps (see 
supplementary information) the result of the complemen- 
tarity protocol 2 — Id-\-i and the smallest eigenvalue of the 
partially transposed state Min [eig(p^^)] . This allows 
us to derive realistic error estimates. In Fig. |4] we plot 
explicitly the matrix elements of the bound entangled 
state {qi = —0.07, q2 = —1.73, qs = —0.5774) together 
with the theoretical prediction. The numerical values 
are 2 - J^+i = -0.032 ± 0.003 and Min [eig(p^^)] = 
+0.0101 ± 0.0017. This clearly proves that the state is 
bound entangled: non-separable, but PPT. 

We studied also the historical Horodecki-state,!^ and we 
can indeed experimentally confirm bound entanglement 
thereof for a large range of the A-parameter. For instance, 
for A = 3.5, we get 2 - i^+i = -0.025 ± 0.003, and 
Min [eig(p^^)] = +0.012 ± 0.002. 

Maximum Complementarity Protocol: In a next 
step we want to verify entanglement of the state, 
shown in Fig. [4j directly by using the complementarity 
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protocol. For that, the correlation functions had to be 
measured directly, which is considered to be difficult 
due to normalization: In quantum state tomography, 
proper normalization is imposed during numerical 
search; here we have to obtain probabilities directly from 
experimental coincidence rates. We obtain the proba- 
bilities P from the normalized coincidence rates F via 

PAk,Bk {h, ik) = ^Ak,Bk {h, ik)/ ^sk,tk ^Ak,Bk {Sk.tk). 

From this we calculate the correlation functions C and 
obtain: 



Correlation function 


Theory 


Experiment 




0.675 


0.667 ±0.005 


Ca2,B2 
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0.463 ± 0.005 


^A3,B3 


0.468 
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Ca4,B4 
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2 - (^4 = ^k^Ak,Bk) 


-0.079 
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The uncertainties are determined from multiple exper- 
imental runs. Experiment and theory agree well, the 
maximum complementary protocol confirms that we de- 
tect a truly entangled state. In literature, criteria for 
entanglement are usually applied on density matrices 
(as we did in the first step). We have performed ex- 
tensive experimental tests and conclude that the max- 
imum complementarity protocol can directly be mea- 
sured if the experiment is well under control. This im- 
plies a significant advance in the exploration of high- 
dimensional entangled photons, since the experimental 
and computational complexity is strongly reduced: We 
need Nqst = dP — ± 4(i^ measurements for (over- 
complete) state tomography, and only N^mcp = d ^ dP 
measurements to determine the function Id^ri directly 
(without normalization), and N^mcp = d'^ -\- d^ includ- 
ing normalization. 

Generality of the Maximum Complementarity 
Protocol: Let us discuss what we would expect for the 
case of higher dimensional 0AM entangled photons in 
the state pd- We search the parameter space of prime 
and prime-power dimensional entangled states (the di- 
mensions where d-\-l MUBs are known), and optimize our 
maximum complementarity protocol to find the states 
where the detection of entanglement via 2 — Id-\-i is 
strongest: 



min 2 - [pd] 
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(8) 



In each case, a similar geometry as shown in Fig. |3] is 
obtained. In dimension d = 2 • 3 = 6 so far only 3 MUBs 
are found and there are strong numerical hints^^ that 
there exist not more; applying only these does not lead 
to the detection of bound entanglement. In Ref.l^ we 
have shown via geometry considerations that if a Her- 
mitian operator detects entanglement of states in a cer- 
tain quantum space W, then it also detects entangle- 
ment in the multi-partite product space W^'^ := {pf^ = 
^Ck,iPk,i\ck,i > O.^Ck^i = 1}. Also in this case, the d'^ 
Bell- type vertex states Pk^i are obtained by applying a 
Weyl operator in one subsystem to Po,o = ^^Pf?- 
Therefore, the maximum complementarity protocol is 
also applicable for multipartite states and can detect 
bound entanglement therein. 

To illustrate the richness and difference between bipar- 
tite and multipartite entanglement, let us discuss first 
the famous case of d = 2. Then the vertex states are the 
Smolin states^^ that are known to be multiparticle un- 
lockable bound entangled, because no local party can dis- 
til entanglement; however, two parties can distil (unlock) 
the entanglement, which would in this case be actually a 
pure maximally entangled state; but it is only available 
for the other parties, not themselves. This case was re- 
cently experimentally demonstrated.^ ^ Surprisingly, en- 
tangled states inside the magic simplex can be dis- 
tilled to the Smolin-type vertex states, which them- 
selves are undistillable. For d > 2, in addition, the mul- 
tipartite magic simplex contains PPT-entangled states, 
therefore, states that cannot even be purified to a ver- 
tex state. They are bound entangled and only "bound 
unlockable" due to their multiparticle nature. 

We have proven experimentally that bipartite bound 
entanglement exists in Nature, which was predicted in 
199^ and started an intensive theoretical quest in un- 
derstanding its meaning. We have chosen orbital an- 
gular momentum entanglement of photons that is scal- 
able in the dimensionality,^^ and our high-quality results 
suggest that further exploration of high dimensional en- 
tangled quantum states is possible. Therefore, our ex- 
perimental method and the maximum complementarity 
protocol is a useful addition to the toolbox to explore 
different types of entanglement including bound entan- 
glement. With the case of two qutrits as the most simple 
system we lay the foundation to pursue the question why 
Nature should provide us with such a strange form of 
highly mixed entanglement that can not be purified, de- 
spite the fact we have used pure maximally entangled 
states as a resource to produce the state: What is its 
physical meaning? What are the quantum information 
theoretic applications of high-dimensional bound entan- 
glement? 

Finally, via the maximum complementarity protocol we 
provide an alternative proof that the maximum number 
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of mutually unbiased bases (MUBs) cannot be more than 
d-\-l: If for families of states that are optimally detected 
by /d+i, another MUB would be added, we would detect 
separable states as entangled ones. The minimum num- 
ber of existing MUBs is known to be the smallest of the 
prime factors of d-\-l. It is also open whether one always 
needs all MUBs to detect bound entanglement. Con- 
sequently, investigating inseparability problems opens a 
different trail to look for a solution of the number prob- 
lem of MUBs. 

Methods: 

Maximum complementarity protocol. As an explicit ex- 
ample we show our choice of basis vectors of the four 
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To detect (bound) entanglement, we choose the following 
combination of correlation functions = |^)*): 

2 

Ca,,b, = ^Tr(|n,mod(n + 1, 3)*)(n, mod(n + 1,3)*^ 

2 

CA„B,=Y.^r{\ik,il){ik,il\p3), k = 2,3,4- 

2 = 

Experiment. We generate the spatially entangled pho- 
ton pairs by collinear Type-I SPDC in a PPKTP crys- 
tal (length L = 2 mm) of a LGq laser beam (Kr+, 
A = 413 nm, beam waist at crystal Wp = 325 ^m, 80 mW 
power). The temperature is tuned to detect a similar 
amount of downconverted photons in the £ = and 
^ = ±1 0AM mode. We image the crystal surface 
with 7.5 X magnification using a telescope onto the SLM 
(Hamamatsu X10468-07) surface. The SLM is operated 
under an incident angle of 10 or 5 degrees; this allows 
us to use a single SLM for both signal and idler photon. 



The (phase-corrected) SLM is used with a blaze angle of 
1 mrad. The elimination of the zeroth order is done in a 
time-reversed fashion before the SLM with a pinhole in 
the far field of the crystal behind LI (Fig. |2|. The far 
field of the SLM surface is sent to the single mode fibre 
using 10 X objectives, with a detection- fibre mode waist 
at the SLM of 1275 ^m. The fibres are connected to 
single-photon counters and we post-select photon pairs 
by coincidence detection (time window 7 ns). All mea- 
surements were integrated for 2 s; we obtain typically 
30,000 single counts and 1,500 coincidence counts for 
conjugated-field detector settings and <10 coincidence 
counts for the cases where no coincidences are expected. 
The SLM kinoforms required for 2-state superpositions 
involve only phase modulation, however, for 3-state su- 
perpositions such as those used for implementation of the 
maximum complementarity protocol, require spatial am- 
plitude modulation, too. Because this is very inefficient if 
done in a continuous way,^^ we employ binary amplitude 
modulation by removing the blaze at spatial positions 
where the normalized modulus of the detection field am- 
plitude is smaller ^/O^. For tomographic reconstruction, 
we obtain an overcomplete set of measurements by de- 
tecting all two-state superpositions with relative phases 
of {0,7r/2,7r,37r/2}. 

Supplementary information: 

Phase space of two qutrits. Below, in Fig. [5) the finite 
phase space of the magic simplex^^ is drawn for dimen- 
sion d = 3. This space is spanned by all nine Bell states 
Pfe^i, which can be generated by applying the Weyl oper- 
ators Wk,i onto one arbitrary maximally entangled Bell 

state, denoted by Po,o: Pk,i = Wk,i 1 ^0,0^^2,/ ^ 1- 
Due to the group structure of the Weyl operators certain 
mixtures of the Bell states Pk^i are geometrically equiv- 
alent, i.e. have the same properties concerning separa- 
bility, bound entanglement, free entanglement and nonlo- 
cality (violation of a given Bell inequality) . Thus one has 
not to analyze all possible mixtures since some are equiv- 
) alent concerning the properties we are interested in. In 
particular, the so called lines are special. A line is formed 
by choosing one Bell state, e.g. Po,o^ and applying to it 
d — 1 times the same Weyl operator, e.g. Wq^i. Apply- 
ing another time the Weyl operator brings one back to 
to original state (in our case P^o)? because of the peri- 
odicity of the Weyl operators. In our case the three Bell 
states {Po,0 7 ^o,i7 ^0,2} form a line in the phase space 
(see Fig. [5|. We find d -\- 1 lines with different orienta- 
tions through one Bell state Po,o (red lines in Fig.|5|. To 
each (red) line there are 3 parallel lines. In summary, 
we find (3 + l)x3 = 12 lines, which have the same ge- 
ometry regarding separability, bound entanglement, free 
entanglement, and nonlocality. 

For our state ps (Eq. 7), we see that the states weighted 
by qi and q2 form a line and qs is the coefficient for an- 
other line that is parallel to the first one. Thus, we have 



7 



li • • • • 




(00) (01) (02) 



• • • • 



FIG. 5. Illustration of the finite discrete classical phase space 
for dimension d = 3 of the locally maximally mixed states of 
the magic simplex W. Each point (kl) represent one of the 
nine Bell states Pk,i- All possible complete lines through the 
point (00) for d = 3 are drawn, representing one class of 
states which have the same geometry concerning separability, 
bound entanglement, free entanglement and nonlocality, i.e., 
those that are unitary equivalent. The same holds for each 
line which is parallel to any red line. 



(3 + 1) X 3 possibilities to choose a line and 3 possibilities 
to weight them with qi and ^2- And there are 2 possibil- 
ities to choose parallel lines to the chosen one weighted 
by qs. Thus we have in total (3 + 1) x 3 x 3 x 2 = 72 uni- 
tary equivalent possibilities to obtain the state ps, and 
all of these have in theory the same geometry regarding 
separability, bound entanglement, free entanglement and 
nonlocality. 

In the experiment, the nine Bell states are not fully uni- 
tary equivalent due to wave front errors, thus we get for 
all 72 possibilities slightly different values that allows 
to perform statistical analysis to obtain the true experi- 
mental errors for the complementarity protocol 2 — I^-^-i 
and the minimum eigenvalue of the partial transpose 
Min [eig(p^^)]. 
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